function[u0]=pdeic(x)

%u0=[1.3*abs(sin(pi*x))+0.6*abs(sin(5*pi*x)) ; 0.5 * abs(sin(3*pi*x))+1.1*abs(sin(5*pi*x))];
u0=[1.3*abs(sin(pi*x))+0.6*abs(sin(pi*x)) ; 0.5 * abs(sin(pi*x))+1.1*abs(sin(pi*x))];
%u0=[1;0];
%u0=[0;1];
%u0=[sin(pi*x);sin(pi*x)];

